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Abstract 

We develop a symplectic method of quantization of lightcone QCD. We find that 
boundary gauge fields are crucial for a consistent and complete quantization. By applying 
the symplectic Faddeev-Jackiw method, we very carefully remove unphysical degrees of 
freedom and obtain the true phase space and the complete Hamiltonian. The result is 
important for the high energy QCD evolution and for a further extension of the JIMWLK 
equation, for which we find the second-order correction. Finally, we make an important 
note about a peculiarity of four space-time dimensions. In additional, our method sheds 
new light on the lightcone quantization of a scalar field. 
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1 Introduction 

This paper is devoted to studying of a high energy scattering amphtude in QCD from 
ab initio. An important ingredient of high energy QCD is the Jahhan-Marian-Iancu- 
McLerran-Weigert-Leonidov-Kovner (JIMWLK) equation [H [2]. Now, there exists the 
problem to make a step beyond the first-order calculations on which the JIMWLK equa- 
tion is based. The primary aim of the paper is to obtain the complete structure of 
quantum Hamiltonian at lightcone, which allows to study the all orders. Recently, its 
importance was emphasized in Ref. [3] where it was used for the high-order studying of 
the high energy evolution in QCD. Although in Ref. [3] the right quantum Hamiltonian 
was used, the rigorous proof is not known. The long-standing problem of quantization is 
not solved yet. In this paper we solve this problem and show explicitly how to obtain the 
complete Hamiltonian and the Hilbert space. 

To quantize the theory we must identify the degrees of freedom. In lightcone theories, 
it is the well-known fact [4] that boundary conditions at x~ = oo are crucial for physical 
consequences. Namely, choosing a different boundary condition, we obtain a different 
physic. Hence, in addition to the bulk fields, we must consider the boundary degrees of 
freedom. 

There exist many fundamental works concerning the lightcone field theory [5l|6]. The 
problem of zero modes is also studied well [7]. However, we assert here that a complete 
solution is not given yet. In a scalar theory, there is no natural boundary condition 
because such theories are too poor. The main idea of Ref. [7] is to enumerate the whole 
set of boundary conditions and to take only conditions giving a consistent Poisson algebra 
with suitable quantization rules. This way has two big problems. 

The first big problem is that there is no natural choice of a boundary condition. 
Since a boundary condition is crucial for physical implications [1], the theory becomes 
physically incomplete. It is known that, studying a lightcone theory, one must treat each 
type of boundary condition separately [8]. Thus, a boundary condition must be intrinsic 
property of a lightcone theory, not being chosen by hand. In the QCD case, since the 
field Ai is a gauge potential, we cannot simple set it to zero at the x^-infinity. Hence, 
the important task is to find a physical assumption that allows us to impose a boundary 
condition on the field Ai at the infinity. The typical physical assumption is to select 
field configurations with finite energy [HE]- In the massive scalar case this easily gives 
the constraint (p{oo) = 0. So, boundary fields do no affect on the physical sector. In 
the QCD, the situation is more complex. It was shown in Ref. [9] that asymptotic of a 
gauge field contains a nontrivial and important information about the structure of QCD. 
Unfortunately, in Refs. [21 UO] only the antisymmetric boundary condition is used. In this 
paper we prove that this is a wrong assumption. Although this choice gives a consistent 
quantization, but it misses an important physical information about the theory. In this 
paper we perform a more careful analysis and show how correctly to obtain boundary 
constraints and to quantize the theory. 

The problem of boundary conditions is closely related with the problem of zero modes 
which causes the appreciable theoretical interest O |7]. In is known that zero modes 
emerge only if the symmetric boundary condition is adopted [7]. In this paper we show 
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that there are no zero modes in the lightcone QCD, since neither symmetric nor antisym- 
metric boundary conditions are apphcable to the real QCD. Fields at the infinity must 
be considered as independent variables. This means that the method used in Ref. [5] 
is inapplicable to the real world. Instead, we show that nontrivial aspects of the theory 
are the complicated structure of the phase space and the Hamiltonian containing the 
infinite number of terms. So, our result disagrees with the Hamiltonian calculated in Ref. 
[5] due to the different boundary condition. Note that the results of Ref. |5] is widely 
exploited in many works. As an example in the context of the high energy scattering, 
higher-order calculation was performed in Ref. [11] where running coupling corrections 
to the JIMWLK equation was calculated. 

The second big problem is the inability of the conventional canonical method to treat 
fields at boundaries. This problem is well known, the classical definition of the Poisson 
bracket fails, since the naive variational derivative does not exist. In Refs. |121ll3j it was 
proved that it is possible to generalize the Poisson bracket to a case with boundaries. 
Although these generalizations seem challenging, currently there are no successful appli- 
cations to physical problems. The situation becomes even more difficult if a theory has 
constraints. Second class constraints are handled by a construction of the Dirac bracket 
which can be properly defined only in a finite-dimensional case. In an infinite-dimensional 
case, due to the existence of surface integrals, there exists the problem to define the vari- 
ational derivative. The non-existence of the variational derivative was used in Ref. [H [7] 
to conclude what a boundary condition gives an inconsistent theory. However, this is 
not the case. Instead of the canonical method, in this paper we employ the symplectic 
Faddeev-Jackiw method [14J which allows to handle fields at the infinity. While naive 
variational derivative can be consistently defined only with zero boundary conditions, 
the symplectic method does not use it. A symplectic structure have no problems with 
the appearance of surface integrals. Applying the symplectic method to QCD and not 
imposing any boundary condition, we remove unphysical degrees of freedom and obtain 
the phase space with the complete Hamiltonian. 

The master logic of this paper is: (1) To establish the physical importance of the gauge 
field at the ^"-infinities. (2) To declare the necessity to quantize the theory without any 
boundary condition. (3) To conclude that the traditional canonical method fails to do 
this. (4) To develop the symplectic method that can quantize the theory. (5) Using 
the symplectic method, to quantize the theory. (6) To check the consistence by the 
rederivation of the JIKWLK equation. 

The key result of the paper is the complete Hamiltonian of QCD. The phase space 
of the theory is the space of fields Ai{x~,x) that obey Ai{+oo,x) = —Ai{—co,x). The 
symplectic form is uj = J d-dAf A dAf. The boundary field 7f (x) is determined by the 
two following equations: 

= pabclH + / {gfabcd^A'^A'^ " 5^+) dx' 
— oo 

where is an external current. The original gauge fields Ai is related to Ai as 

Ai{x~,x) = Ai{x^,x) - ^7i(x) 
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The momentum (x , x) is determined by condition vr^ (±00, x) = and the constraint 

a_7r- + d.d.At - gfabcd-A\Ai + 5J+ = 
The Hamiltonian is given by 

H[A,] = i(vr-)2 + \f^.[A,]F:^[A,] - gAUl 

This paper is organized as follows. In Sec. El using the standard canonical method, 
we briefly review basic properties of the lightcone QCD. In Sec. [3l we formulate the 
symplectic Faddeev-Jackiw method. In Sec. [H we use this method to quantize a scalar 
field with the boundary condition which will be useful in the QCD case. In Sec. [5l 
we apply the symplectic method to the lightcone QCD. This section is divided on nine 
subsections: (|5.ip the primary analysis of the phase space, ()5.2p the first step of the 
Faddeev-Jackiw algorithm and the complete set of Gauss' constraints, (j5.3p the analysis of 
the gauge invariance, (j5.4p the gauge fixing, (j5.5p the residual gauge transformations and 
the finite-energy condition, (j5.6p the formulation of the complete Hamiltonian, (j5.7p the 
quantization of the theory and the construction of the Hilbert space, (15. 8p the perturbative 
expansion of the Hamiltonian operator, (j5.9p the third-order expansion which shows how 
the theory lives in space-times having more than four dimensions. In Sec. [6] we use the 
first-order Hamiltonian to rederive the JIMWLK equation. Then we obtain the 0{a1)- 
correction to this equation. Section [7] contains our conclusions. 



2 QCD at lightcone. 



This section contains the standard canonical approach to QCD at lightcone. Also, we 
recall typical quantities and fix the notation. 
The covariant derivative is defined as 

D^^l^ = {d^ - igA^)i^ (1) 

where if) \s a. color multiplet. The gauge curvature tensor is defined from the commutator 
of covariant derivatives 

ig-\D^,D,]^j = F^,il, (2) 

F^u = d^A^ - dyA^ - ig[A^,A^] (3) 
In the subsequent calculations we shall sometimes use the the adjoint index notation 

\ = A^T, (4) 

[Ta-,Tb] = ifabcTc (5) 

Sp{Tan) = (6) 

C = Qy^K - QvAl + gfab.A^'^Al (7) 
The Lagrangian density of QCD is 

L = --Fl^F^^ ^ gAlJ>^ (8) 
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where Jj^ is a current of matter fields. For the fermion case it has a form 

Jj^ = ^I^^Tai; (9) 

The equations of motion is 

5L _ dL dL 



Straightforward calculations give 

dL 



d{d,A-) 



-F7 (11) 



8T 

- = gfabcFrK + 9J^ (12) 



^\ = + gfabcFrAl + gJ^ (13) 



A shorthand notation is 

i| = D.Fr+5j„^ (14) 

In the lightcone coordinates, the canonical momenta are 

dL 



Ft (15) 



There is a primary first class constraint 

< = (16) 

Since F*"*" do not contain the velocities, the standard lightcone second class constraints 
are 

<-i^r=o (17) 

where i = 1,2. The Hamiltonian density is 

H = 7r^d+A''^-L (18) 

The fundamental Poisson brackets are 

{A'^^ix),7ray)} = 6'j^^6{x-y) (19) 

Since the Hamiltonian (|18p is defined only on the constraint surface which is the image of 
the Legendre map, an extension of the Hamiltonian on the whole phase space should be 
given. Usually, such extended Hamiltonian is called total Hamiltonian. The consistency 
condition requires 

Ut,H} = (20) 

This leads to 

= ^ = D,F:^+ + gj+ = Z?,< + gj+ (21) 



where we have used the definitions (jl4p and (|15p . Eq. (j2ip is the new first class constraint 
which is usually known as Gauss' law. 

Now, we impose the first gauge fixing condition 

= (22) 

This converts the constraint ([160 to a second class constraint. In this gauge we have 

= d+A^ (23) 

H = F''+F+^ + i (2F+^F+^ + 2F-iF-' + FjF^^) - gA^.Ji' = 

^F^^+F+^ + ^F.iF-' + iF.jF'^ - gA^J^^ = (24) 
\ {F-+F+^ + F'+F+i) + \F^^F-' + - gA^J^' = 

Using the fact 

F'+ = -Fi_ (25) 

we obtain the final result 

H = ^{7T-f + ^F,,F^^ -gA.Jf^ (26) 

where we have dropped the obvious color summation index a. 
The second gauge fixing condition is 

= (27) 

This converts the constraint (j2ip to a second class constraint and the momenta 7r~ be- 
comes a depended variable. The constraint (fT7|l has been converted to 



vr* - d-Ai = (28) 

where we have used the identity F'"*" = F^i and the gauge fixing ()27p . Similarly, in this 
gauge the constraint ([2T]) can be expressed as 



a_7r- + - gfabcO-A^A^ + 5J+ = (29) 



The next task is to express vr" from (|29p and to substitute it into the Hamiltonian (I26p . 
Finally, one can obtain the Hamiltonian as a function of only the fields Af with the 
canonical structure derived from the constraints ()28p via the standard Dirac brackets. 

It is instructive to obtain the free version of the Hamiltonian by setting the coupling 
constant to zero: g = 0. Ignoring boundary terms, we have (tt")^ = (diAi)'^ and this 
term is canceled by the equivalent term within FijF^K So, we have 

Ho = \diA]diA] (30) 

This Hamiltonian can be easily recognized as a sum of ordinary massless scalar field 
theories in the lightcone, which Hamiltonian is FL^ = ^(Sjc^)^. Hence, to quantize the 
free QCD, we can take the standard result of quantization of the free scalar theory. The 
index j in ()30p directly corresponds to transverse polarization of a physical gluon. 
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3 Faddeev-Jackiw method 



In this section we briefly review a simplified and geometrized version of the symplectic 
Faddeev-Jackiw method [14J . The method is a useful tool to study constrained dynamical 
systems. 

Consider a mechanical system whose Lagrangian is 

L{qiAi) = aMW -y{(l) (31) 

Since it is linear over the velocities , we have a constrained system. This system can be 
studied by both the standard canonical method and the original Faddeev-Jackiw method 
[14j . However, if one needs observables that depend only on the coordinates q*, we can 
directly endow the coordinate manifold Q with a symplectic form. A Lagrangian L{qi, qi) 
can be viewed as a function on a tangent bundle. In each fiber, the expression ai{q)q^ 
of the Lagrangian ()3ip defines a linear form. Hence, on the whole manifold we have the 
diff'erential 1-form ^ = ai{q)dq^ . The exterior differential d produces 2-form 

^ = di = ^dqi ^dq' = -(^-^^ dqi A dq' (32) 
^ dqi ^ ^2 \dqi dq' J ^ ^ ^ ^ 

where A is the wedge product. The form lo is closed: dco = due to nilpotency = 0. 
To be a symplectic from, a closed 2-from must be nondegenerate. If w is nondegenerate, 
the phase space is completely specified and a Poisson bracket is generated by the inverse 
of the skew-symmetric matrix coij 

The Hamiltonian is simple 

H = V{q) (34) 

If oj is degenerate we must make the additional steps. A degenerate 2-form to at each 
point of the manifold defines its own kernel in the linear tangent space. For any vector 
a G Keruj and for any tangent vector 6, we have uj{a,b) = 0. It is easy to check that 
the Lie bracket of two vector fields a,b £ Keroj is also an element of Kercj. This fact 
is a consequence of = and the Cartan formula for duj{a,b,c) where a,b G Kerw. 
Hence, due to Frobenius theorem, Kerw is an integrable set of vectors and it at each point 
defines a submanifold. So, if there are no topological obstructions, we have a foliation 
of the original coordinate manifold which becomes a union of parallel submanifolds of 
smaller dimension. A fiber of the foliation will be also called Kerw. The foliation gives 
an equivalence relation: two points are equivalent to each other if they belong to a 
same submanifold. The main statement is: there exists a manifold Qi and a natural 
isomorphism i such that 

QixKera; = Q (35) 

Also, the isomorphism i induces a closed and nondegenerate 2-form uji on the manifold 
Qi, since Kerco is factored out. Hence, the form oji is a symplectic form on Qi. For 
practical purposes, instead of a usage of the formal statement (j35]) it is sufficient to find 
coordinates that explicitly separate Kercj from Q. Let {zi} be coordinates in Kerw and 
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{(jk} be the remaining coordinates. In this coordinates the form co does not depend on 
the forms dzi 

LO = idki{q)df A dq^ (36) 

Also, it should be stressed that ujkiio) cannot depend on the variables due to the 
property du = 0. So, the variables form the manifold Qi endowed with the symplectic 
form u)i = u}{q). Note that in the symplectic analysis the decomposition ([35]) and the 
well-known Darboux's theorem have the same origin. 

Since we have Qi = Q/Kerco, the manifold Qi is a quotient manifold, while the 
Hamiltonian H is a function on Q. Hence, in general cases H cannot be reduced on Qi. 
This is the most intimate point of the story. To obtain a reduced Hamiltonian, we can 
try to fix a natural left inverse of the canonical epimorphism Q — t- Qi. The fixation is 
based on the fiber diffeomorphisms whose vector fields lie in Kercj. Namely, we construct 
a condition 

dH{KeT uj) = H{q, z + Sz) - H{q, z) = (37) 
In the coordinate notation, we are looking solutions of the equation 

Taking coordinates q^ of a point in Qi, we solve Eq. (j38p and find the remaining one z*. 
The reduced Hamiltonian will be Hi(q) = H{q, z{q)). Nevertheless, we have problems 
again. Since, in general, Eq. (I38p can have: no solutions, precisely one solution, more 
than one solution, and any z* as a solution. Since the last alternative arises if H does 
not depend on z*, we can restrict H on Qi trivially. If there are no solutions, we must 
remove such points q'' from the theory. Hence, we have new constraints. More precisely, 
if Eq. (I38p has a solution only in a submanifold Q2 C Qi, then Q2 contains the true phase 
space of the theory. Currently, we are not concerning the case of more than one solution, 
since we shall see that in the application to QCD only the first and last alternatives are 
required. 

Finally, after the solution of Eq. (I38p we have a submanifold Q2 with Hamiltonian 
H2 = Hi. Again, if Q2 7^ Qi we have a problem. The reduced on Q2 symplectic form wi 
may be degenerate. We denote it uj2. If UJ2 is a degenerate 2-form in Q2, then, to remove 
Kera;2, we just repeat the procedure described above. This gives new pairs (Qsif^s)) 
{Q4, UJ4) and so on, until at some step n we shall find the true phase space that is a triple: 
a coordinate manifold Qn, a nondegenerate closed 2-form Wn, and a Hamiltonian Hn that 
is a function on Qn- Now, the method is completely specified. 

It is quite remarkable that the presented method can be safely applied to systems 
having infinite and continuous number of degrees of freedom, because functional deriva- 
tives are not required. There is only one potentially problem to invert a symplectic form, 
but it can arise only at the last step of the construction. 

The method is also useful to impose an external constraint. Since a differential form 
can be naturally restricted on a submanifold and the Hamiltonian is 0-form, we are able 
to impose by hand any constraint which defines a submanifold in the original manifold. 
A dynamical system will be completely specified if a reduced symplectic form is non- 
degenerate. How to consider the degenerate case it has been described above in this 
section. External constraints are required to impose a gauge fixing condition. Since there 
is freedom to choose a gauge, in Sect. [5] we shall use a manually constructed external 
constraint. 
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4 New lightcone quantization of a scalar field 



In order to understand the role of fields at the boundary we consider a scalar field theory 
in 1 + 1 dimensions. The free Lagrangian is 

L = y d-ipd+Lpdx~ (39) 

Since L is linear over the velocities we are able to apply simplified Faddeev-Jackiw 
method which was constructed in Sec. [3l Since the coordinate space is a linear space, 
the tangent space is naturally isomorphic to it. The Lagrangian L induces 1-form 

^ = [ d^ip{x)dip{x)dx (40) 



where dip{x) is a functional analog of 1-forms. It acts on a vector u[x) as d(p{x) (u) = u{x). 
Also, L induces a symplectic 2-form u in the linear space of all fields (f 

u = d^ = j d-d(p{x) A d(p{x)dx (41) 

where A is the continuous wedge product. Taking any two fields ipi{x~) and (p2{x~), the 
form uj can be expressed as 

U){ipi,Lp2) = / (5(/7i(/72 -<fid(f2)dx~ (42) 



where we have used a shorthand notation d = d-. The constant symplectic form u 
promotes the space of field to a phase space. Note that if there are no restrictions on 
the field values at the infinity, the form u is not degenerated because one cannot find a 
field (fii such that the action of the form (j42p is zero at any ip2. Only if we impose the 
symmetric boundary condition (^(— oo) = ip{+(X)), then w degenerates in the space of field 
configuration that is symmetric at infinity because uj(const, ip) = for any ip. In this case 
there exists a so-called zero mode - the constant field which commute with each other 
field. This is the typical result of Discretized LightCone Quantization (DLCQ) with 
the symmetric boundary condition. Nevertheless, in general, the symmetric boundary 
condition can have physical meaning only in a nontrivial topology case. 

With aim to apply the result to QCD we impose the following boundary condition 

ip{-oo) = (43) 

and there are no restrictions on (p(+oo). The condition (I43p will arise in QCD due to the 
opportunity to impose a gauge fixing condition. 

To obtain a Poisson bracket we have to find an inverse of the matrix w. However, 
a direct inversion meets with math difficulties. A careful inversion requires some math 
stuff which will be developed now. 

Let us define three linear spaces. Let Q be the space of all fields if{x~). Let Qo be 
the space of fields obeying the constraint (f43l) . Let Q be the space of fields obeying the 
antisymmetric boundary condition (^(— oo) + <p{+oo) = 0. There is an important fact 
that Qo and Q are linear subspaces of Q. This means that these subspaces are invariant 
under any linear combination aipi + b(p2- As an example, constraints like (p{—oo) = 1 
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do not generate linear subspaces. The symplectic form oj can be naturally restricted to 
a linear subspace. In general, a restricted form can be degenerated, but in our case it is 
easy to check that the corresponding forms in Qq and Q are not degenerated. Hence, the 
spaces Qo and Q have a natural phase space structure. The key step of our setup is a 
construction of a linear map X '■ Qo ^ Q 

x(¥'o) = ^0- ^V'oC+oo) (44) 

where 999 ^ Qo- It is clear that Ker^ = and Imx = Q. Moreover, there exists the 
unique inverse map '■ Q ^ Qo 

X-^{^) = ^ + ^{+<x) (45) 

where (p ^ Q. Hence, we have proved that the map x is an isomorphism of vector spaces. 
This fact allows us to use x as a change of variables in the dynamic system. Hence, we 
shall use the space Q as new coordinates of the system. Also, the map x is a symplectic 
map (symplectomorphysm) . This means that for any (p,^ € Qq we have 

^{xM, XW) - ^{^, V') = ^'/'(oo) J dip- ^V'(oo) J dip = (46) 

So, in the sense of Hamiltonian mechanics, X is a canonical transformation. 

At first sight, the map ()44p seems strange. To see it as a linear map we replace the 
whole x~-line [—00, +00] by a finite interval [—L, +L] and approximate by a finite lattice 
(pi = ip{xi), i = . . . n, xq = —L, Xn = +L. The finite-dimensional linear space formed by 
the set {(fi} is an approximation to the space Q. The map (jH]) becomes (pi = ipi — lpq/2. 
In the matrix notation we have (pj = {6ji — 5oi/2)ipi. Thus, we see that the map (I44p can 
be approximated by a linear matrix Cji = 6ji — 6oi/2. 

Having the space Q and the symplectic form uj (|42p . we are ready to find a Poisson 
structure. This can be achieved by the search of a linear operator P{x — y) such that 

(^(j^i, J P{x - y)ip2{y)dy^ = j (pi(p2dx (47) 

for (pi,ip2 G Q- For brevity, we just show here that a standard ansatz P{x) = as{x) 
motivated by the well-known results of lightcone field theories solves the task. Here we 
have used a sign function 



1 


if 


X 


> 





if 


X 


= 


-1 


if 


X 


< 



e(x) = < if X = (48) 
[ -1 

Its derivative is e'{x) = 25{x). Using the rule 

o +00 

d(p{x)e{x — y)dx = I d(p{x)dx — I d(p{x)dx = —2(p{y) (49) 



we find the solution of Eq. ()17|) 



P{x-y) = -^-e{x-y) (50) 
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The final step is a quantization of the theory. This is the standard textbook result so 
we write only the final result for completeness. Restoring the transverse coordinates, we 
declare the quantization rules 



[ip{x ,x),(p{y ,y)] = --e{x -y )5{x - y) (51) 
^(x-,x)= / -^{ale^^^^--^'^^ + a,e-^^^'^-^^^')^^ (52) 

fe+>0 

= 5{k~ - p~)5{k - p) (53) 

where the restriction > is imposed since the free Hamiltonian must be bounded 
from below. 

A price for using the map (j44p is that the Hamiltonian should be recalculated in 
terms of fields (p. Using the inverse map (|45p . we have H = H{(p + (p{+oo)). This is a 
potentially dangerous expression, since we must define a meaning of the operator (^(+oo) 
and to check commutation relations between observables and the Hamiltonian. In Ref. 
[3] a similar operator is called zero mode. However, ip{+co) is not a zero mode. It is just 
a tool to calculate right expressions for observables such as a Hamiltonian, preserving 
friendly commutation relations. 



5 Complete Hamiltonian of QCD at lightcone 

In this section we apply the Faddeev-Jackiw method [14J to QCD. 

5.1 Bare phase space 

Bare phase space is the space that one obtains after the first step of the Faddeev-Jackiw 
algorithm. Separating the velocities from the QCD Lagrangian density ([8]) which at 
lightcone has the form 

L = \f+J + F+iF_i - ^F.jF'^ + gA^,P (54) 

we identify the following schematic structure 

L = Li + L2d+Ai + L:id+A^+'^{d+A^f (55) 

where Lj do not depend on velocities. There is no dependence on 9+^+, there is a linear 
dependence on and there is a typical quadratic dependence on dj^A^. So, the 

simplified Faddeev-Jackiw method cannot be directly applied. To obtain a phase space 
and a symplectic structure we have to perform the Legendre map only over the variable 
d+A_ 

H = TT-d+A- - L (56) 

where the terms, having the velocities d+Ai, have been removed into the canonical 1-form 
by the way which was shown in Sec. El and the velocities d-^-A- are extracted from the 
identity 

vr- = F+_ (57) 
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So, the Hamiltonian density is 

H = \{^a? + (d-K - gfabcAXAt) + (d^AX - gfabcA\Af) F-^ + -^F-rJ - gAlJ^ 

(58) 

The hnear phase space is the set of fields 

{A+,A_,7r-,^J (59) 

The canonical 1-form is 

i = J n-dAt + J F^.dAt (60) 

where the integration over the space coordinates x~, x' is assumed. The canonical 2- form 
is 

uj = dC = J diT- AdAt+ J dF^^ A dAf (61) 
Note that cj is a degenerate form since, at least, it does not contain dA^. 

5.2 The complete set of Gauss' constraints 

We have seen that the 2-form (j6T]) does not depend on dA^ and A^. Hence, the choice 
of coordinates (j59p explicitly separates Kerw and Eq. (j6ip already has the form of Eq. 
([36|) . So, the coordinate A^ plays the role of variables of Sec. [3l To obtain the reduced 
Hamiltonian, we have to impose a functional extension of the condition (j37p . Making a 
variation 5A^, we have 

= J {H[A+ + 6A+] - H[A+])d^x (62) 

For fixed A-,tt~ , Ai, Eq. (I62p has no solutions or any A^ as a solution. Hence, in 
accordance with the discussion in Sec. [3l we have a constraint within the bare phase 
space. The bulk part of the variation (I62p gives the ordinary Gauss' law (12ip 



D.TT- + D.F^i + gj+ = (63) 

while the boundary one gives two new constraints 

vr" [x~ = ±oo,x) = (64) 

F_i (x~,|x| =oo) =0 (65) 
and the reduced Hamiltonian density is 

H = liTT-f + \Ft^Fi^ - gAtJ- - gAUi (66) 



The constraints (j64p and (|65p are the necessary and important part of the theory. In 
particular, we shall see below that the constraint ()64p is very helpful in the construction 
of the true phase space. 

The constraints (I63p and ()64p completely determine the coordinate vr". It will be 
useful to find vr" and to substitute it into the Hamiltonian ([66]) . However, one should be 
careful, since we shall see that the constraints ([63l) . (|64]l also restrict the coordinates Ai. 
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5.3 Gauge invar iance 



Until now, we did not use any gauge fixing procedure. In principle, we can avoid this 
step by the introduction of unphysical degrees of freedom into the scene. Nevertheless, 
when using noncovariant methods, a gauge fixing allows to simplify a calculation. How 
to motivate a gauge fixing within the context of the symplectic method? In the canonical 
approach, the Dirac procedure says that in a presence of first class constraints we can 
impose new gauge fixing constraints that convert all constraint to the second class. In 
the symplectic method a gauge invariance can be handled in a more straightforward way. 
Let us recall the definition of a gauge transformation. It is a matrix-valued function 
U{x~,x^) £ SU{N). In the matrix notation the gauge field transforms as 

= UA^U^ - -{d^U)U^ (67) 
9 

The Hamiltonian (|58|) is invariant under U by the intrinsic gauge invariance of the theory. 
The Gauss' constraints (|63p. ()64p . and (j65p is explicitly invariant too, since they transform 
covariantly. 

Since we have used the differential calculus, it is instructive to review a some geometric 
stuff. The primary fields j4_, j4j} at fixed 2;+ forms an infinite-dimensional manifold. 
The gauge transformation ()67p acts like a diffeomorphism. Since a diffeomorphism can 
be naturally extended on the tangent bundle, any gauge transformation has the following 
action on a tangent vector: 

= UA^U^ (68) 

where the function A^ implements a tangent vector. Eq. (j68p is a consequence of the 
well-known fact that a difference between two connections is an adjoint vector under a 
gauge transformation. Also, the transformation (j68p naturally induces the action of the 
gauge group on a cotangent vector. Namely, the differentials dA transforms as 

dA^ = uUA^U (69) 

Eq. (169p is a consequence of the requirement that the convolution dA^{A^) must be a 
gauge invariant. Similarly, since vr" is an element of the cotangent bundle, it transforms 
in according to (j69p . So, we can conclude that the 2- form (I6ip is a gauge invariant. 



5.4 Main gauge fixing 

Since the Hamiltonian, the Gauss' constraints, all physical observables, and the canonical 
form cj are gauge invariants, the gauge invariance (j67p is a genuine invariance of the 
dynamical system. So, we have the opportunity to impose a constraint on the set of 
fields (j59p within the bare space Q. Mathematically, this is just a choice of a submanifold 
Qi C Q. Certainly, we must prove that such a choice is physically admissible. 
We use the standard lightcone gauge 

= (70) 

It is easy to show that each field configuration has a gauge equivalent obeying (j70|) and 
there are no Gribov ambiguities. Hence, the choice (I70p does not restrict the physical 
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sector of the theory. A gauge transformation that leads to the condition (|70p is the usual 
parallel transporter 



U{x ,x) = Pexp(—ig J A^{x ,x)dx j (71) 



where the coordinate Xq of the initial plane is arbitrary. 
After the gauge fixing, the canonical 2-form (I61|) becomes 



00 = J dF^i AdA'} = J d-dAI A dAI (72) 



which is just a vector version of the form ()4ip in a scalar theory. The Gauss' constraint 
([63|) becomes 

9_7r- + d.diAl - gfabcd-AlAi + gj^ = (73) 

Now, after the gauge fixing and the elimination of vr^, our phase space is entirely 
specified by the set of the transverse gauge potentials Ai obeying the constraint (f65]) . 
This space is endowed with the closed 2-form ()72p and the Hamiltonian 



H = \{^-a? + \Ft,F:^ -gA'iJi (74) 



We mentioned above that the Gauss' constraints also impose additional constraints on 
the variables Ai. To see the hidden set of constraints, we integrate the Gauss' constraint 
r3|) over the space variable x~ . Using Eq. we have 

diAi{+QX), x) = diAi{-<x, x) + G[Ai, J] (75) 



where J] is the functional containing the bulk integration. The constraint (j65p 

becomes Ai{x~ ,oo) = 0. Again, integrating over x~ we have a restriction on field 
values at the infinity 

j4j(— oo, oo) = Aj(+oo, oo) (76) 

Thus, we see that ^j(iboo,x) are not independent variables. More precisely, only a field 
Ai obeying (f75]l . ([65]) . and (|76l) is an element of the phase space. 

In principle, at the current stage, one might finish the determination of the phase 
space. However, to construct a quantum theory, one should invert the symplectic form in 
the nonlinear manifold formed by the field Ai obeying the Gauss' constraints. Without 
additional assumptions, this is a difficult task. Typically, the naive restriction ^i(oo) = 
is used or the more refined antisymmetric boundary condition Aj(+oo) = — Aj(— oo) is 
adopted O [10]. However, we have seen that all such conditions are physically irrele- 
vant. To obtain a tractable quantum theory having a complete physical information, we 
continue the further detailing of the phase space. 



5.5 Residual gauge fixing 

Although the residual gauge freedom is a well-know fact |10] . it cannot be rigorously 
studied within the canonical method, where, at first, one has to identify first class con- 
straints. Thus, the canonical method fails to reproduce a true phase space, since one has 
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to use value of a field at the infinity as a canonical variable. In the symplectic method, 
we can safely use a constraint that involves field at a boundary. 

The gauge fixing (j70p is not complete. There exist gauge transformations that does 
not change In fact, any U{x), which does not depend on x~, is a symmetry of the 
dynamical system. So, after imposing the constraint (j70p . we again have the opportunity 
to impose an external constraint that does not miss physical information. To obtain 
a gain from a residual gauge fixing, we have to decide a constraint that will be used. 
One might try to use the residual gauge freedom to make the antisymmetric boundary 
condition Ai{—oo,x) = —Ai(+oo,x). However, in general, this is impossible by the 
following reason. Consider the transverse plane with two gauge connections Ai(x) and 
Bi{x). We are trying to find a gauge transformation U{x) that leads to Ai = —Bi. Since 
U{x) must exist for any Ai[x) and Bi{x), it is sufficient to find a counterexample. If 
Bi = Ai, then Ai = Bi for any U. Hence, it is impossible to make an antisymmetric pair, 
until the case Fij [A] = is in question. Since a connection with zero curvature can be set 
to zero by gauge transformation, in this case the simple counterexample does not work. 
The case of Fij = will be processed later. 

It is easily to show that the residual gauge freedom cannot make the symmetric 
boundary conditions, too. Indeed, the condition Ai = Bi and Eq. ()67p immediately give 
the restriction Ai = Bi which is obviously violated if Ai and Bi are arbitrary fields. 

There is one more source of external constraints. We can restrict ourself to use 
only field configurations whose energy is finite. More precisely, we are free to impose 
constraints that remove only field configurations having infinite energy. It is clear that, if 
Fij(zizoo, x) 7^ 0, the Hamiltonian energy (I74p is infinite. So, we impose the new constraint 

F,j-(±oo,f) = (77) 

This means that at x~ = ±00 the field Ai{x) must be a pure gauge field. Although 
such fields have the trivial holonomy, it is again impossible to make the antisymmetric 
boundary condition from the residual gauge freedom. There exists an intrinsic alge- 
braical obstruction. Let Ai and Bi be gauge fields. To make the antisymmetric boundary 
condition, one has to find a gauge transformation U that leads to Ai + Bi = 0, where 
transformed fields Ai and Bi was given in Eq. ()67p . In the matrix notation, we have the 
equation 

U{Ai + B,)U^ -2-{d^U)U^ = (78) 

g 

It has a solution only if the field {A-\- B)/2 is a pure gauge field. Hence, a solution exists 
only if 

j,.,\A + B- 



2 (79) 

Using the definition ^ and the constraints Fij [A] = Fij [B\ = 0, we can recast Eq. (j79P 
to 

[A,Aj] + [B,,Bj] - [Bi, Aj] - [A, Bj] = (80) 

To construct a counterexample, we choose Bi = 0. In this choice, Eq. (I80p becomes 

[A,,Aj]=^ (81) 

If a theory has a non-Abelian gauge group and if it has two or more transverse dimensions, 
then there exists a pure gauge field that violates Eq. ([8T]) . So, the inapplicability of the 
antisymmetric boundary condition has been proven. 
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Instead, we impose the following realizable constraint: 

^i(-oo,f) = (82) 

It is indeed realizable. Let Ai(x) = Ai(—oo,x) be a field such that = 0. Let xq 

be any point in the transverse plane. Let U{x,xo) be a gauge link from the point xq 
to the point x. Since the field Ai has zero curvature, U{x,xo) does not depend on form 
of a path, it only depends on the initial and final points of the paths. To achieve the 
constraint (I82p . we perform the gauge transformation V{x) 

V{x) = U-^{x,xo) = U{xo,x) (83) 

It is easy to check that after the gauge transformation (|83|) every gauge link becomes 
trivial. Also, it is important that there is no the Wu-Yang ambiguity [15j which says that 
two distinct gauge potentials can have the same field strength. The strong constraint 
Fij = completely removes such complications. So, we have proven that the gauge 
choice (j82p is allowable and it completely fixes the residual gauge freedom. 



5.6 The Hamiltonian 

Now we are ready to formulate the main result. At first, let us collect the obtained facts. 
Let 7j(x) be value of a field Ai{x~ ,x) at the boundary x~ = +oo 

ji{x) = Ai{+oo,x) (84) 

Due to Eq. (f77|) and ([75]) . 7^ is not independent variable. It obeys the equations 

+00 

da-{x) = J {gfahcd-A\Ai - gj+) drc" (85) 
—00 

- djlf + gfabclhj = (86) 
Moreover, Eqs. (j65l) and (I82p give the asymptotic of and Ai at the transverse infinity 

^i(x-,oo) = (87) 

Note that the term fabcd-A^A'^ in (I85p is just the gluon contribution to x^ component of 
the overall color current. Thus, the constraint (j85p has the form dijf{x) = gp^{x) where 
p"'{x) is the transverse density of color charge. 

Now, we apply the method which was presented in Sec. HI After the choice of the 
coordinates ^ 

Ai{x' ,x) = Ai{x^ ,x) - --ii (f ) (88) 

we have 

Ai{+oo,x) = -Ai{-oo,x) (89) 

i,(+oo,x) = ^7i(f) (90) 

+00 

= \fabclH + / {gfabcd-AlAi - gj+) dx" (91) 
—00 
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H = \{^a? + \Ft,F^J -gA^l (92) 
9_7r- + d.diAl - gfabcd-A\Ai + gj+ = (93) 
uj = j d-d A'} A dA'} (94) 

We shall treat Eq. (|9Up as a constraint that is defined in the space of pairs {Ai,ji). In 
this space, the constraints (j86p and (|9ip completely fix the field 7^ in terms of the field 
Ai. Note, that the Hamiltonian (|92p is uniquely defined only on the submanifold that is 
formed by the three constraints ([86]) . (f90|) . and ([9T]) . 



5.7 Quantization 

Unfortunately, the attempt to directly quantize the theory fails due to the lack of a 
set of coordinates that forms a linear space and solves the constraints. In Sec. S] we 
saw that if a phase space, having a continuous number of degrees of freedom, can be 
parameterized by a linear space, than the symplectic procedure can be safely applied 
(a similar parametrization is just the isomorphism x in Eq. ()44p ). So, to quantize the 
theory we have two distinct ways: either to solve the constraint or to introduce unphysical 
degrees of freedom. We adopt the second way. 

If we remove the constraint (I90p . the phase space can be linearly parameterized by the 
set of fields Ai, since the remaining constraints fix only the field ji. The symplectic form 
(|94p is already defined in the proper subspace. After the fixation of 7^, the Hamiltonian 
H becomes a complicated functional over the field Ai. Moreover, one can add to H 
any functional that vanishes on the constraint surface (j90p . The symplectic form a; is a 
nondegenerate form both with the constraint ()90|) and without it. 

Now, the quantization is straightforward. The phase space is just tensor products of 
the phase space of the scalar field theory for each transverse index i and each color adjoint 
index a. Applying the rules ([5TH53]) . we can construct the Hilbert space in the basis of 
the free theory which is the limit g = 0. Also, in the free theory, there are no unphysical 
degrees of freedom, since it is easy to check that 7^ = and the constraint (j90p gives 
Ai(oo,x) = 0. Dropping the constraint (|90p . we extend the phase space from the zero 
boundary condition to the antisymmetric boundary condition. At lightcone, it is believed 
that these two boundary conditions give the equivalent quantum theory and one can safely 
use the extended phase space. Indeed, for the both cases the commutation relation ([STIl 
does not change. However, in general, the quantum equivalence is not proven. Although 
it is well known that, in comparison with the free case, the complete Hamiltonian with 
an interaction can significantly change the structure of the Hilbert space, we hope that 
for practical purposes of high energy evolution in QCD and the perturbation theory it is 
sufficient to work in the free Hilbert space where possible unphysical degrees of freedom 
do not affect observables. 

As we have seen, the symplectic method works well for linear spaces. In a more gen- 
eral case, to quantize a theory, a self-consistent Poisson algebra should constructed. In 
addition, if a second class constraint exists, then the Dirac bracket should be found to 
obtain new constraintless Poisson algebra. Indeed, from Poisson viewpoint the constraint 
([90|) is a second class constraint. Finally, the algebra should be quantized and appropri- 
ate Hilbert space with a Hamiltonian bounded from below should be constructed. The 
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problem is complicated by the existence of functionals involving spatial derivatives and 
nontrivial boundary conditions. In this paper, we just conjecture that the constraint (|90|) 
can be safely removed. 

5.8 Perturbative expansion 

The Hamiltonian (j92p has the very complicated structure as a functional over the funda- 
mental fields Ai. Moreover, it contains arbitrary large power of the coupling constant g. 
Indeed, 7j, having finite order of g, cannot be a solution of Eq. ([55]) . For purposes of the 
perturbation theory it is useful to expand H in the series on g. Here we limit ourself to 
the first-order calculations and assume that the typical parton number has order 1 (in a 
dense case it can have order g~^). 

In the zero-order ji = and we just obtain the free Hamiltonian (j30p . Using Eq. 
([88|) . in the first-order the constraints ([86]) and (j9T]) give 

da- - dj^f = (95) 
= J {gfabcd-A\Ai - gj+) dx- = gp-{x) (96) 

— oo 

Eq. ()95p says that the 1-from jidx^' is closed. Since the transverse plane with the boundary 
conditions (|87p has trivial de Rham 1-cohomology, any closed 1-form is an exact form. 
This means that there exists 0-form such that 

7i = diip (97) 

Hence, the solution of Eq. (f96]l is 

= gdr^P^i^) = IJ ^^,P%y)dy (98) 

where the integral representation has been written for two transverse dimensions. In the 
first order Eq. ([93]) becomes 

5_7r- + d.diA'i + gjt = (99) 

where we have defined the density of the overall color current 

j+ (x- , x) = -fabcd-A\Ai + J+ (100) 

which obviously obeys / j^(x, x~)dx~ = Pa{x). To solve Eq. ([99]) . we define the operator 
dZ^ by the following way: 

j^fi^') = 11 - y')f{y')dy- (loi) 

~oo 

Other choices of dZ^ are also allowed, but all they will give the same result. Recalling 
the constraint (jM]) and the boundary condition ([89]) . we are able to write the action of 
dZ^ on Eq. ([99]) 

t:- + d,A'i + g^3t = ^ (102) 
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Substituting vr^ from (jl02p into the Hamiltonian (j92p and again using Eq. (j86p . in the 
first-order we have 

= i^o + + 5i-if (103) 

where jj is the transverse components of the overall color current 

f ) = -f,,,d^A';A'^ + Jt (104) 

Note that in the derivation of the Hamiltonian (I103P we did not use the solution (I98p 
for 7j. In fact, in the first-order it is not required. Nevertheless, it will be required in 
the expansion up to the second-order. Second-order calculations is more complicated and 
lies beyond the scope of this paper. For this case we only briefly sketch the calculations. 
Let be the second-order correction to 7^. Fortunately, before the third-order, Eqs. 
(I95|),(I96]), and ^ is still vahd also for ^j. Eq. ([99]) is also vahd, but instead of Eq. (fTOUl) 
we have 

j+(x-, x) = -fabcd-AlAl - fabcd-A\^t + J+ (105) 

Finally, one can substitute vr" from (jl02p into the Hamiltonian (|92p and expand the result 
up to the second-order of the coupling constant g. However, the explicit form is tedious 
and it is preferable to filter simultaneously the terms in according to a matrix element in 
question. 



5.9 Third order 

The expansion up to the third order is much more complicated, but it allows us to see the 
announced peculiarity of four space-time dimensions. Let r/j be the third-order correction 
to 7j. The constraint (f86l) gives 

dtri^j - djrit = -gfabdhj (106) 

Now, the 1-form rj = rjidx^ is not closed. Using the language of differential form, we can 
rewrite Eq. (jl06p in the very compact form 

dVa = <Ja (107) 

where aa is the given 2-form: aa = —gfabcJiljdx^ A dx^ . Eq. (|107p immediately gives 
the integrability condition 

daa = (108) 

It is remarkable that in two transverse dimensions the condition (jl08p is fulfilled auto- 
matically! If there are more than two transverse dimensions, then this condition must be 
checked explicitly. If it is not fulfilled, then Eq. ()106p does not have a solution and the 
whole theory fails to be consistent! More precisely, in this case the finite-energy condition 
()77p cannot be imposed and the theory should be constructed in a different manner. We 
shall not touch this problem and just assume daa = 0. 

In the general case, the solution of ()106p can be formulated with the help of the 
Poincare duality. There exists the Hodge map * that maps from n-forms to {D — n)- 
forms, where D is the dimension of the transverse plane. The Hodge map defines the 
Hodge adjoint of the exterior derivative: 5 ~ *d*. Since the 1-cohomology is zero, the 
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harmonic 1-forms do not exist and the Hodge theorem states that any 1-form is a sum of 
exact and coexact 1-forms 

rj = dtp + 6d (109) 

where ip is a 0-form and ^ is a 2-form. The exact part dip disappears from (|107p as weh 
as in the first order. So, Eq. (I107P becomes 

d5e = a (110) 

If dcTa = then a solution for 9 obviously exists. The coexact form 59 disappears from 
Eq. ([96]) . since it is easily to prove the identity di{59)i. So, diip can be directly found from 
the extension of Eq. (|98p where now p"'{x) is a more complex function over Ai,J~^ ^j. 

If the transverse plane is two-dimensional, the Hodge decomposition (I109p has the 
more simple form 

r]i = ditp + Ejidj-ip (111) 

where 'iIj{x) and ^(x) are two arbitrary functions. The tensor Sij is the alternating tensor 
such that Eij = —Eji and £12 = 1. The constraints ([861) and ([911) uniquely determine 
and via the following equations: 

Alp = Fi (112) 
AtP = F2 (113) 
where A = didi is the Laplace operator and Fi, F2 is functions depending on Ai,J~^,ji, 

The fourth and higher orders can we processed in the way that has been developed 
for the third order. 



6 JIMWLK at once 

In this section, to check our result, we straightforwardly derive the JIMWLK equation 
[2] from the first-order Hamiltonian. The general method of derivation was developed 
in Ref. |16|. It is argued that only one soft gluon emission into new opened phase 
space contributes to the evolution of the elastic scattering amplitude, since only a linear 
term over small 6Y survives in the differential equation. Although it is proved that 
one needs to known boosted wave function of a projectile with second order accuracy of 
the Hamiltonian perturbation theory, the key building block is the soft gluon emission 
amplitude of the first-order. 

In Ref. |16j the different gauge condition is used. To use the Hamiltonian (jlOSp . we 
have to recalculate the matrix element of the interaction part. Only the term diA^dZ ja 
gives a contribution to this matrix element. The operator diA^ emits a soft gluon, while 
the operator probes the fast color charges in the projectile. In the first-order, the 
wave function of the boosted projectile with a soft gluon can be calculated with the help 
of Eq. (21) in [l6]: 

|p(i)) = -i7o-iy|p(°)) (114) 

where V = J diAfdZ^ja- Using the definition (jlOip and the gauge version of (j52p . we 
need to calculate the integral 

e^'=+^"ie(x- -y-)(ix- = e^''^^' f e'^^''~le{x~)dx- (115) 



2 ' " ' J 2 
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where is the longitudinal momentum of the soft gluon. The integral in Eq. (jllSp can 
be calculated with the help of the following integral identity for any t > 0: 

+00 

J e^*^^e(x)dx = J (116) 
—00 

Then, we obtain the integral 

"e"'^y~jt{y-,x)dy- (117) 



Since the operator j+ probes the color charges having a high longitudinal momentum, 
in the evaluation of the integral in Eq. ()117p we can neglect the small longitudinal 
momentum of the soft gluon. Hence, the integral (I117p can be expressed as 

/ e'''^y'j+{y-,x)dy- = f j+{y-,x)dy- = pa{x) (118) 
J k+=0 J 

The action of the operator Hq^ is 2k'^ /}? as it is shown in [16j. Collecting all terms we 
obtain from Eq. ()114p the first-order correction to the projectile wave function 

where the integration over k'^ is performed in the new opened rapidity window. It is 
useful to define the gluon creation operator in a given transverse point 



at,+ . = e-'^^al ^ -^=^ (120) 



d^k 

^x,k+,i,a.~ J ^ "'k,k+,i,a ^(27r)2 

To express Eq. (|119p in the term of Eq. (|120p , we have to calculate the following integral: 

^e'^^d'^k = 2m-^ (121) 
k r 

Finally, at the first-order we have 

Eq. (I122|) has a clear physical description. A fast color source at a transverse point x 
can emit the soft gluon into arbitrary point y of the transverse plane. 

The further calculation is fully equivalent to Ref. [16]. If the projectile scatters in the 
target color field aa{x), the evolution equation for the elastic scattering amplitude is 

-d^=(h^ 1^) (123) 

^ ' ^ V +2VUz)p\y)Sp-{x) I 



_ {z-y)iz-x) 
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where Vbaiz) is the gluon scattering amphtude which is a Wilson hne in the adjoint 
representation. The operator S is the operator of S'-matrix. 

Let us briefly consider the general properties of high-order corrections. The detailed 
study will be a subject of a separate paper. It should be emphasized that we assume 
that the number of partons in the projectile has the order 0(1). This differs from the 
approach of Ref. [3] where the order 0{l/g) was assumed. This distinction gives different 
corrections to the JIMWLK equation. Thus, a physical meaning of each approach should 
be clarified. 

The careful examination of the terms in the Hamiltonian gives the following form for 
an operator that has a nonvanishing contribution to the JIMWLK equation: 

Ai Dij^] (125) 



5_ 

where D[j^] is a polynomial over the current . Indeed, the operator Ai creates only 
one soft gluon, while all other operators must act on the fast partons. Only the operator 
gives nonvanishing matrix element in the high-energy limit. The integral kernel dZ^ 
is required to compensate the multiplier in the numerator, which comes from Hq^ in 
the perturbation theory (I114[) . The convolution of the space and color indexes in (11251) 
is not restricted. 

Analyzing the Hamiltonian (192^ , we conclude that the form (1125^ arises only from the 
term (vr~)^. In the first order we already used the corresponding term diAfdZ^ja- In 
the second order there are no terms having the form (I125p ! In the third order only the 
following term gives a contribution to the JIMWLK equation: 

V = -g^AhUabcj-jt (126) 



where 7^ is given by the first-order expression (j98p . The calculation of the wave function 
of the boosted projectile is very similar to the first-order calculation. Starting from Eq. 
(jll4p . we obtain the following correction to the wave function 

I (1). _ [ ig'fabc {Xl-X2y \Z-Xi\\z-X2\ ^ . n'' ( X,) d'' ( Xo) ' IP^) fl27) 

where ^ is an infrared regulator and the integration is over xi,X2, z,k^ . Note that in 
higher orders there exists the problem of operator ordering in the Hamiltonian. We 
adopt here the standard prescription: it is assumed that the Hamiltonian is already 
normal ordered H =: H :. From practical viewpoint, such normal ordering means that 
two operators p°'{xi) and p'^{x2) cannot act on a same parton. Thus, in Eq. (1127^ the 
normal ordering is required. Applying the algorithm of Ref. |16j . we obtain the following 
0(a^)-correction to the JIMWLK equation: 



^ -{5,{p«(y),f«(xi,X2)}} \ 

+2VU^p\y)Sf-{xi,X2) 
\ +2VUz)f\xi,X2)Sp''{y) j 



\P) (128) 



{z-y){xi-X2) , \z-xi\\z-X2\ , . 
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f^(fl,f2) = fabc ■ P''{xi)p''{x2) 



(130) 



where the bracket {, } denotes the anticommutator. Note that Eq. (jl27p does not give the 
0(a^ )-correction to the JIMWLK equation, because, in addition, one must also examine 
the 0(5r^)-terms in the Hamiltonian. 

7 Conclusion 

In this paper, we have constructed the complete lightcone Hamiltonian of QCD. It has 
the very complicated structure and infinite power over the coupling constant g. However, 
the perturbative expansion can be performed in a controlled way. A further task is to 
apply the developed method to study of various physical problem. Specifically, one can 
systematically calculate higher corrections to the JIMWLK equation (sometimes, they 
are called NLO, NNLO and so on). 

Unexpectedly, our calculations give a restriction on the dimension of the space-time. 
If the space-time is three-dimensional, the residual gauge freedom can be directly used 
to make the antisymmetry boundary condition. Indeed, in one transverse dimension Eq. 
(|79p holds automatically. Hence, there is no reason to introduce boundary field ji. So, 
three-dimensional QCD has the much more simple Hamiltonian which contains a finite 
number of terms and a finite power of g. It is natural to expect that this theory is too 
insignificant in comparison with the four-dimensional QCD. 

If the space-time has more than four dimensions, a quantization of the theory should 
be different. Probably, the finite-energy condition (I77p has no solutions. This means that 
all excitations of such a theory have infinite energy. Although now we cannot prove that 
higher-dimensional QCD does not exist, the exceptional role of four-dimensional QCD is 
clearly visible. This result is not very surprising. The history of quantum field theory 
knows enough examples where the theory chooses a geometry. 

For completeness, the two-dimensional QCD should be mentioned. The developed 
method does not give any profit in this case. The reasons are: the field Ai does not exist, 
the residual gauge freedom degenerates to global color rotations, there are no propagating 
degrees of freedom - there are no gluons, and obviously high energy evolution is trivial. 
However, the theory is not fully trivial - it is still possible to study the g^-potential and 
Wilson loops. 

Let us make several remarks about the method of lightcone coordinates. Although 
after the first studies this method rises great expectancies, today it is clear that the 
method is more complicated. While in the Cartesian coordinates the complexity of the 
theory is pushed entirely into the quantum level, in the lightcone coordinates it remains 
in the classical level via the complicated structure of the phase space and the uncom- 
putable Hamiltonian. Due to the "conservation" of the complexity, there are no reasons 
to expect a simple solution in the both approaches. However, in the lightcone method the 
complexity grows in a controlled manner. Indeed, in the classical level we work with the 
classical mathematical background which is much more friendly than the Hilbert space. 
Possibly, since the lightcone theory has the trivial vacuum, the solution of many problems 
can be found in the classical level. Since the Lagrangian can be generalized to a fully 
covariant form, there is no a priori frames to quantize the theory. Hence, the lightcone 
coordinates is a natural choice as well as a choice of the Cartesian coordinates. 
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